Non-adiabadic charge pumping in a hybrid SET transistor 
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We study theoretically current quantization in the charge turnstile based on the hybrid (SINIS 
or NISIN) SET transistor. The quantization accuracy is limited by either Andreev reflection or by 
Cooper pair - electron cotunneling. The rates of these processes are calculated in the "above-the- 
threshold" regime when they compete directly with the lowest-order tunneling. We show that by 
shaping the ac gate voltage driving the turnstile, it should be possible to achieve the metrological 
accuracy of 10~*, while maintaining the absolute value of the quantized current on the order of 30 
pA, just by one turnstile with realistic parameters using aluminium as superconductor. 

PACS numbers: 73.23.Hk,74.45.+c,84.37.+q 



Nanoscale tunneling structures provide the general ba- 
sis for development of metrofogical sources of electrical 
current utilizing controlled transfer of individual charges 
However, despite the beautiful achievements based 
on experiments with gated arrays of metallic tunnel 
junctions [1, H, 0, and with semiconductor surface- 
acoustic- wave and charge-coupled devices [1, 0, B ll] , no 
fully satisfactory system in terms of both the accuracy 
and current magnitude has been realized yet. It was sug- 
gested recently [l^l that an unexpectedly simple struc- 
ture, a single-electron (SET) transistor with two hybrid 
normal metal - superconductor (NIS) or superconduc- 
tor - normal metal (SIN) tunnel junctions holds promise 
as a quantized source of current. The first experiments 
with such a transistor as a turnstile [lo| demonstrated 
correct operation at the level of classical charge dynam- 
ics, but they were not yet conclusive as to its ultimate 
accuracy. In this Letter we analyze theoretically all the 
relevant higher-order quantum tunneling processes which 
limit this accuracy. The main conclusion we reach is 
that these errors can be suppressed in a single ordinary 
aluminium-based device to the level mandated by the 
metrological requirements (< 10~*), while keeping the 
absolute current relatively large (see Fig. 2] below), pro- 
vided the single-electron charging energy of the turnstile 
is sufficiently high. 

The basic "classical" dynamics of the hybrid SET tran- 
sistor (Fig. [T]) as a charge turnstile can be described 
conveniently on the stability diagram shown in Fig. [21 
Periodic variation of the gate-induced charge ng{t) = 
CgVg{t)/e with frequency / (notations are defined by 
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FIG. 1: (a) Hybrid SET transistor with SIN or NIS tunnel 
junctions, and (b) time dependence of the ac gate-induced 
charge rig = CgVg/e oscillating with frequency / around the 
point rig — 1/2. 



Fig. [T|) indicated by the line with arrows in Fig. [5] drives 
the transistor periodically between the two nearest sta- 
bility areas, e.g., n — and n = I, where n is the equi- 
librium number of extra electrons on the island. The 
turnstile operation requires that the lowest-order tunnel- 
ing transitions are organized so that at finite bias volt- 
age V and low temperature T they transfer prec isely 
one electron per period 1/f through the transistor [lO|. 
The properties of the tunneling thresholds (solid lines 
in Fig. [2) that make this possible in the hybrid transis- 
tor but not in the normal-metal one can be seen from 
Fig. [21 The thresholds in the hybrid are shifted with re- 
spect to the normal-metal system (dashed lines in Fig. [2]) 
by the superconducting energy gap A, i.e., the shift along 
the Jig axis is (5 = A/2Ec', where Ec = e^/2CE and 
Cs — Ci + C2 + Cg, expanding the stability areas. As 
a result, the neighboring stability areas overlap, and the 
gate voltage can drive the system between them keep- 
ing it all the time in the region of suppressed tunneling. 
Also, in this process, when the outgoing gate voltage tra- 
jectory leaves the initial stability area, it crosses only one 
of the tunneling thresholds that define this area, allowing 
electron tunneling in only one direction. For instance, if 
the state n = is brought by increase of rig out of its ex- 
panded stability area into the n — 1 area (Fig. [2]), elec- 
tron can tunnel into the transistor island only through 
the left junction. When the gate voltage decreases back 
to n = state, electron can tunnel out only through the 
right junction [l^]. 

This turnstile operation is possible for any, e.g. har- 
monic, periodic time dependence ng{t) with the ampli- 
tude sufficiently large to move the system between the 
two stability areas (Fig. [2]). The time that the system 
spends, however, in the overlap region of the two areas 
does not play any useful role in the turnstile dynamics, 
and on the contrary, increases the effect of the unwanted 
transitions. In order to maximize the turnstile opera- 
tion frequency and the output current, one needs then to 
minimize this time by making the waveform ng(t) as in 
Fig. [TJd. In this case, the system is switched abruptly be- 
tween the regions where electron tunnels in or out of the 
transistor, and the operation frequency / is limited only 
by the need to make the probability of missing these tran- 
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FIG. 2: Charge stability diagram of the hybrid SET transis- 
tor operated as a turnstile. Dashed lines are the tunneling 
thresholds of the rhombic stability regions n — 0, 1 in the 
normal-metal case. Solid lines show the thresholds in the hy- 
brid transistor shifted by the superconducting gap A. Peri- 
odic variation of the gate voltage (line with arrows) transfers 
one electron per period through the transistor. 

sitions e~'^l'^^ sufficiently small. At zero temperature, the 
corresponding tunneling rate is ^(JJ) = 7o(C/^/A^ — 1)^/^, 
where \J is electrostatic energy change due to tunneling, 
7o = G/S./e^ and G is the junction tunnel conductance. 
Optimized waveform (Fig. [TJd) should be abrupt on the 
time scale of the turnstile period 1//. It should, however, 
be smooth on the scale /i/A to avoid excitations of the 
higher-energy states of the transistor leading to errors in 
the turnstile dynamics. This condition can be satisfied 
easily, since for a typical current of 100 pA, the frequency 
/ = //e < 1 GHz is wen below A//i ~ 50 GHz. 

In addition to missed cycles of tunneling due to fi- 
nite frequency /, the basic correct tunneling sequence 
can be interrupted by thermal excitations due to finite 
temperature T, or quantum higher-order tunneling pro- 
cesses [ll| which set the theoretical limit on the accuracy 
of the quantized current I — ef produced by the turn- 
stile. The rate of thermal errors depends on how far 
the gate-voltage trajectory is from the crossing points 
of the four relevant tunneling thresholds shown as solid 
lines in Fig. ^ The thresholds are given by the condi- 
tions Up = A on electrostatic energy change Up due to 
forward (wanted) or backward (unwanted) electron tun- 
neling in the jth junction: 

Ut = ±2Ec{vi+ng-l/2), Up = ±2Eciv2-ng + l/2) , 

where Vj is a part of V that drops across the jth junc- 
tion: vi = (C2 + Gg/2)V/e and V2 = (Ci + Cg/2)V/e. 
These equations show that at the thresholds of cor- 
rect tunneling, the energy barriers for unwanted transi- 
tions through the opposite junction of the transistor are 
A — Ui = A — U2 = eV. Thus, with exponential accu- 
racy, the thermal probability of electron tunneling in or 
out through the wrong junction leading to no net charge 
transfer in the cycle, is g-eV/fesT^ Another type of un- 
wanted thermal transitions is the excitation of an extra 
electron through the transistor during the part of the pe- 
riod spent in the overlap region of the two stability areas. 
Electron is transferred by two successive excitations over 
the energy barriers A — Up , so that the thermal excita- 




FIG. 3: Tunneling thresholds for several different tunneling 
processes: Andreev reflection (AR) and Cooper-pair/electron 
(CPE) cotunneling, in a hybrid SET transistor. As in Fig. [2] 
solid lines are the thresholds of the lowest-order tunneling. 
All tunneling processes are driven by bias voltage V and are 
allowed above the corresponding threshold. 

tion exponent for the overall process is g-i^^-f'^)/i^BT ^ 
Comparing the probabilities of the two types of thermal 
errors, we see that the thermal error rate is minimum for 
eV ~ A: in practise, the resulting classical error e~^/'=sT 
is less than 10~® at realistic temperatures T ~ 100 mK. 

We consider now quantum errors assuming ideal s-wave 
BCS superconductors in the hybrid transistor structure. 
The rates of "elastic" higher-order processes which trans- 
fer electrons coherently, without creating excitations in 
the electrodes, are different in the NISIN and SINIS 
structures. In the NISIN transistor, the dominant elas- 
tic process is electron cotunneling, the rate of which is 
smaller than the rate 7 of the lowest-order tunneling 7 
roughly by a factor {hG/e^){SE/ A) [H, where SE/A is 
the level spacing of the transistor island. For typical pa- 
rameters, e.g., ^m-size island, this suppression factor is 
very small, about 10~^ — 10~^, but does not quite reach 
the metrologically required level. In the SINIS transistor, 
in the relevant regime eV ~ A, the main contribution to 
elastic leakage is due to rectification of the ac Josephson 
current through the transistor. The resulting dc current 
is proportional to the square of the SINIS critical current 
and is much smaller than the inelastic leakage assisted by 
Andreev reflection that is considered below. 

The rates of incoherent "inelastic" processes depend 
only on the local properties of the tunnel junctions and 
are the same in the NISIN and SINIS transistors. Inten- 
sity of these processes decreases rapidly with the num- 
ber of involved electron transfers. The simplest process 
of electron inelastic cotunneling through the transistor 
is energetically forbidden in the relevant voltage range 
eV < 2 A. Transitions next in the order of complexity are 
Andreev reflection (AR), i.e. tunneling of two electrons 
in a Cooper-pair for which the superconducting gap does 
not provide an energy barrier, and Cooper-pair/electron 
(CPE) cotunneling. Electrostatic energy gains in these 
processes are 

AR: U++ = 'iEcivi+ng-l), U++ = 4Ec{v2 - Ug) , 
CPE: Wp = Up + eV, (1) 
and the diagram of the corresponding thresholds, U — 
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for AR, and = A for CPE is shown in Fig. [3l If single- 
electron charging energy is small, Ec < A, (i.e. 6 > 1/2 
in Fig. [3]) AR is allowed in the regions of the lowest- 
order tunneling needed for turnstile operation. Each AR 
process causes an error by transferring one uncontrolled 
extra electron. For larger charging energy, Ec > A, turn- 
stile can be operated in the regime with suppressed AR 
{6 < 1/2 in Fig. and only the higher-order CPE pro- 
cesses cause errors. Qualitatively, in the CPE, instead 
of one electron jumping in or out of the transistor is- 
land, this transition is combined coherently with electron 
transfer of another electron through the whole transis- 
tor. To avoid creating superconducting excitations, the 
necessary tunneling of two electrons in one of the transis- 
tor junctions in this process happens as AR. This CPE 
cotunneling is allowed energetically for any turnstile pa- 
rameters and limits the accuracy of current quantization. 

Quantitatively, we calculate the rates of the two 
higher-order tunneling processes assuming the simple 
quasi-lD ballistic geometry of the turnstile junctions, in 
which different transport modes in the electrodes are not 
mixed by tunneling. This assumption is reasonable in 
view of large conductivity of electrodes of practical SET 
transistors. Because of the non-adiabatic variation of the 
gate voltage (Fig. [T]d), both higher-order tunneling pro- 
cesses take place in the "above-the-threshold" regime, 
when they coexist with the lowest-order single-particle 
tunneling. We start with the rate jab. of the Andreev 
reflection. Above the single-particle threshold, the stan- 
dard description of AR as the two-step transition pertur- 
bative in the electron tunneling amplitudes t (see, e.g., 
[isj) should be modified to account for the competing 
single-particle tunneling with rate "f{U~^). Similarly to 
the theory of the Coulomb- blockade threshold , this 
can be done simply by taking into account the lifetime 
broadening ij{U^)/2 of the initial state. 

Because of the mutual coherency of Cooper pairs in 
different orbital states in the superconducing electrode, 
the amplitudes of the Cooper-pair tunneling from dif- 
ferent states p within each transport mode into the two 
single-particle states with energies e^, e; in the normal 
electrode should be summed coherently. The total AR 
amplitude A is then: 



Aiek,ei) 
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+ ek — u 
u = U+ +i-i{U+)/2, 



(2) 



where Up,Vp = [(1 ± ep/flp)/2]^^^ are the usual BCS 
quasiparticle factors and Op — (A^ -I- SpY^^ is the quasi- 
particle energy. Taking the sum over p under the stan- 
dard approximation of constant density of states p and 
tunnel amplitudes t in the relevant energy range on the 
order of energy gap A, we get 

A{ek,ei) = pt'^A[a{u - et) + a{u - e;)]. 



a(e) = (e2-A2)-i/2 



A-e+(e2- A2)i/2 



(g2 „ A2)l/2 



(3) 



The main qualitative feature of the amplitude A is the 
resonance at the gap edge, e ~ A, where the rate |Ap 
diverges as l/|e — A|. Level broadening, in our case due 
to the single-particle tunneling with rate 7, broadens the 
resonance and suppresses the divergence. 

The amplitude A gives the total rate of AR at small 
temperatures ksT <^ A : 
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where in the adopted quasi-lD model the states fc, I in 
the sum should belong to the same transport mode. The 
result of summation over these modes can be expressed in 
terms of the normal-state conductance G within the nat- 
ural junction model in which transparency varies ex- 
ponentially with energy on the scale cq ^ A. The effec- 
tive number J\f of the transport modes in the junction is 
determined then by the decrease of transparency with in- 
creasing transverse energy of the mode: Af — Smeo/Trh"^, 
where S is the junction area and m is electron mass. The 
sum over modes and integration over the total energy can 
then be done separately giving the AR rate: 



JAR 



I&ttN 



defie-U++/2)f{-e-U++/2) 



\J2a{±e + Ec-il/2)\^, g = hG/e^ 



(4) 



If AR transitions are not energetically allowed, the 
leakage current is determined by the third-order CPE 
cotunneling which combines AR with one more electron 
transfer in the opposite junction. The part of the CPE 
amplitude A that corresponds to the two-electron AR 
transfer process is calculated as above for direct AR. 
Combining terms with different ordering of the three in- 
volved electron transfers we get the total CPE amplitude 
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'2Ec + 2U+ -Sk-ei ek + ei - U+ - u' 
Hek - U+) + aiei - [/+)] + [a(u - e^) + a(u - e/)]- (5) 

1 1 



U+ + u - Ek ~ ei 
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U+ -u + eu + ei 



Summing all transitions with this amplitude as above, we 
obtain the total rate of the CPE cotunneling: 



ICPE 



dn- 



deu / dei\A\^ 



[1 - f{^k)][l - f{ei)]f{n + ek + ei-U+- eV). (6) 

Figure shows the gate dependence of the zero- 
temperature normalized rates of the (wanted) single- 
particle tunneling, 7 = 7/70, of the AR transitions, 
lAR = 7Ai?,/(7o3/167rA/'), and CPE cotunnehng, ^cpe = 
7cPB/(7o5^/327r^A/') at the optimum bias point eV — A 
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FIG. 4: Performance of the SINIS hybrid turnstile, (a) 
Normahzed rates (for normahzation, see text) of the vari- 
ous processes as functions of gate position at the optimum 
bias point eV = A: CPE cotunneling (curves in the top 
part), single-particle (lower right corner) and Andreev tun- 
neling (taller curves on the right). The different sets of curves 
refer to -Ec/A = 1 (black dash-dotted), 2 (red dashed), and 4 
(blue solid). The thresholds at i?c/A = 4 for single-particle 
and AR are indicated by vertical arrows, (b) The maximum 
pumped current ([7| as a function of the allowed error rate p 
for i?c/A = 10, 4 and 2 from top to bottom. 



for a few values of the ratio Ecj As in Fig. [3l the 
thresholds of single-particle and Andreev processes coin- 
cide for Ec/^ = 1, but for larger values of this ratio 
there is a window between the two onsets. Kinks in CPE 
rate occur at these thresholds, marked by dashed vertical 
lines in Fig. 0^ for Ec/^ = 4; in between, jcpe changes 
only little. One can see from this plot that an optimum 
gate value - fast single-particle transfer and errors only 
by CPE - exists for the case Ec > A, and it lies within 
1/2 + (2A - eV)/4:Ec < Ug < 1 - eV/4Ec, closer to the 
upper end of this range. The turnstile should thus be 
operated by a gate voltage (Fig. [TJj) switching between 
n„ 



To make a quantitative estimate of the performance 
of the turnstile, we note that since the CPE contributes 
one extra transferred electron, the relative transfer er- 
ror is p = in the operation window discussed 
above. This gives the (maximum) junction conductance 
which can still suppress the CPE error to below p as 
g = 47r[A/'p7/7cP£;]^^^- On the other hand, one can drive 
the turnstile at a frequency / = 7/(21n(l/p)) to suppress 
the missing cycle errors to the same level. The maximum 
current of the turnstile at the error rate p is then 



'MAX 



eA 27r 



1/2 



(7) 



such an Ug and 1 



Og. 



Figure [3)3 shows /max versus p for the most common 
hybrid system using aluminium as the superconductor, 
for which A/fc^ ~ 2.5 K. In this plot, we also take into 
account that J\f oc Eq^ because of the junction area de- 
pendence of both of these quantities, and use an estimate 
of the tunnel barrier characteristics Af — 10^ for Ec = 1 
K. We can see from Fig. [4] that increasing Ecj ^ indeed 
improves the turnstile performance, and a single turnstile 
with Eel ^ = 4 reaches an accuracy of 10~* at about 30 
pA current with ~ 400 kil junction resistance. With 
i?c/A — 10 (such high Ec'-S. were obtained, e.g., in jlSj]), 
100 pA current can be reached with the same accuracy. 

In summary, we have shown that a simple hybrid SINIS 
turnstile should qualify as a metrological source of cur- 
rent. In order to reach sufficient level of current, either a 
very large charging energy or a few parallel turnstiles are 
needed. The latter option is affordable because of the 
simplicity of the basic device In practical pumps 

p^ . other sources of fluctuations (e.g., variations of the 
background charge) that can not be precisely predicted 
by theory, influence the performance as well, although 
the simplicity of our turnstile should make it stable also 
with respect to these fluctuations. 
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